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Abstract. We introduce and discuss a new class of (multivalued analytic) transcendental func- 
tions which still share with algebraic functions the property that the number of their isolated 
zeros can be explicitly counted. On the other hand, this class is sufficiently rich to include all 
periods (integral of rational forms over algebraic cycles). 



1. Motivations. The question about the number and location of (isolated) zeros for 
different classes of analytic functions appears in connection with numerous problems, 
both in geometry (enumeration problems) and analysis. For instance, the famous Hilbert's 
Sixteenth problem (one of the two in the Hilbert list which still remains unsolved) on limit 
cycles of planar polynomial vector fields can be reformulated as a question on isolated 
zeros of the Poincare displacement function (see below for more details). 

The obvious class of functions for which the counting problem admits a universally 
known, simple and precise answer, is that of polynomials (in one complex variable). The 
Fundamental Theorem of Algebra states that the number of isolated complex roots of a 
polynomial of degree d > is always between 1 and d, and is exactly equal to d if the 
roots are counted with multiplicity. 

An instant generalization of this example is given by the class of algebraic functions 
of one variable: the required bounds are provided by the Bezout theorem. Because of the 
ramification and multivaluedness, the counting problem requires certain precautions to 
state. In this particular case it is sufficient to triangulate the complement CP 1 \ S to 
the finite ramification locus 2J, count the number of isolated zeros of each branch of the 



2010 Mathematics Subject Classification: Primary 34C08, 32S65, 34M03, 34M56; Secondary 
34C07, 34M35, 34M99. 

Key words and phrases: Regular flat connections, root counting, Fuchsian systems, monodromy, 

quasiunipotence, differential fields, Picard-Vessiot extensions. 

The paper is in final form and no version of it will be published elsewhere. 



[1] 



2 



GAL BINYAMINI E.A. 



function in each simplex (triangle) of the triangulation, and add up the separate results. 
As with the polynomials, if all roots are counted properly, the exact answer for the total 
number of isolated zeros is given by the intersection theory in purely algebraic terms. 

However, essentially broader classes of such type are not so easily produced. For 
instance, as was argued in |Y05] . no nontrivial class admitting hnite count of zeros, can 
exist as a subclass of single- valued functions on CP 1 \ E with any finite or infinite locus 
E. Indeed, if all excluded points are poles of finite order (and their number is finite), 
then such a function / is necessarily rational. If, on the contrary, one of the points 
a £ £ is an essential singularity for / (the absolute value |/(t)| grows to infinity faster 
than polynomially as t — >■ a), then by the Picard theorem, / assumes almost all values 
infinitely many times in any neighborhood of a. 

Thus any nontrivial class of functions allowing for an explicit global counting of their 
zeros on CP 1 should necessarily include multivalued functions. This necessarily brings 
about the question on the automorphy (monodromy) of these functions, that is, how 
different branches obtained by analytic continuation avoiding the locus E are related to 
each other. 

The most natural assumption, which covers a huge number of possible applications, 
is the assumption that the linear span generated by all branches of these functions, is a 
well-defined finite-dimensional linear space near any point outside E . In particular, this 
includes all algebraic functions whose branches are simply permuted by the monodromy 
transformations. By the classical Riemann-type arguments |IY08j . these spaces can be 
identified with the spaces of solutions of linear ordinary differential equations (or systems 
of such equations) with single- valued coefficients holomorphic on CP 1 \ E. 

If a singular point a £ E is an isolated essential singularity for the coefficients of the 
equation, then one can again argue that solutions of the equation will grow faster then 
polynomially (when approaching a along suitable narrow sectors, to avoid the problems 
with multivalued continuation) . Then by the general principles of the Nevanlinna theory, 
these solutions will have infinite numbers of isolated zeros accumulating to a. Thus with- 
out loss of generality, we may assume that all coefficients of the linear equations have 
only finite order poles, i.e., are rational functions of the independent variable t with the 
polar locus E C CP 1 . However, even this class should be further reduced to expect finite 
count for isolated roots, as explained in the next section. 

2. Counting roots of solutions of linear systems. 

2.1. Crash course on linear systems. We recall a few basic facts on systems of linear 
ordinary differential equations; all these facts can be found in IYQ8, Chapter III]. 

A system of linear differential equations with rational coefficients on CP 1 can be 
defined using a (n x n)-matrix O = {0.^} of rational I-forms on CP 1 , by the Pfaffian 
equations 

n 

dXi — ilijXj, i = I, . . . , n, or in the matrix form dX — flX = 0. (I) 
i=i 

In an affine chart t € C C CP 1 the matrix I-form f2 can be written as A(t) dt with a 
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rational matrix function A(t) — {Ay(i)} so that the Pfaffian system (fTJ) takes the form 
of a system of linear first order differential equations — — = A(t)X. The ("fundamental") 
solution of this system is a nondegenerate multivalued analytic matrix function X = X(t), 
detX(t) for t ^ E, ramified over the polar locus E of the coefficients; the point at 
infinity may or may not be singular. 

The fundamental solution of the linear system ([l} is defined uniquely modulo a con- 
stant right matrix factor X(t) X(t)C, C £ GL(n,C), hence the linear space spanned 
by all entries X{j (t) of any solution, is independent of the choice of this solution. There- 
fore for any simply connected subset T C CP 1 \ E, the linear span is a well-defined 
n 2 -dimensional subspace S£t in the space &(T) of holomorphic functions on T. 

For a closed path (loop) 7 : [0,1] — > CP 1 \ E, 7(0) = 7(1), the result of analytic 
continuation of a nondegenerate matrix solution X(t) along this path is another solution 
X'(t) — X(t)M where M is a constant matrix factor. This factor is called the monodromy 
matrix associated with the loop 7. It is defined by the free homotopy classes of paths 
up to the conjugacy M H> C~ 1 MC independently of the solution X(t). In particular, if 
a, G E and 7 is a simple loop contractible in (CP 1 \ E) U {c^}, then the corresponding 
monodromy is called (modulo the conjugacy) the monodromy of the small loop around 
di G E. 

Example 2.1 (principal). Consider the Euler system with only two singularities, t = 
and t = 00, both of which are simple poles of the coefficients. In the affine chart such 

a system takes the form = — ■ X, where the constant residue matrix A can be 

assumed to have a Jordan normal form (i.e., diagonal). Solution of this system is a 
matrix function t A — exp(^41ni). In the case A is diagonal, the space spanned by the 
components of the matrix X, is the space of linear combinations of the power functions 
Y] j Cj t Xj , j ; = 1, . . . , n, where Ai, . . . , A„ are the eigenvalues of A (in the general case one 
will have also the functions i Aj ln fc t for multiple eigenvalues). 

For a general linear system with rational coefficients there is a well known classifica- 
tion of singular points in terms of the growth rate of the corresponding solutions. This 
classification is parallel to the classification of singularities of single- valued functions into 
poles and essential singularities: if solutions of a system grow at most polynomially near 
a point a £ E, then this point is called a "regular" (moderate) singularity, otherwise the 
singularity is "irregular" (wild). The only caveat requires that the growth estimate for 
multivalued (ramified) functions be valid in sectors with the vertex at a (and not in an 
arbitrary simply connected domain which may be spiraling around a). An elementary 
calculation shows that if a is a simple (first order) pole of the matrix form Q, then it is 
necessarily a moderate singularity, like the points 0, 00 for the Euler system. First order 
poles are called Fuchsian singularities, and the Pfaffian system (|TJ) is called Fuchsian, if 
it has only Fuchsian singularities on CP 1 . In any affine chart a Fuchsian system takes 
the form 





(2) 
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where Ay, . . . , A m are the residue matrices at the singular points ay, ... , a m . The point 
t = oo is singular if and only if Ay + • • • + -<4. m 7^ (in the singular case the residue at 
infinity is the negative of this sum) . 

In general the inverse statement is not true and a moderate singularity may well be 
non-Fuchsian. However, by a suitable local meromorphic gauge transformation a moderate 
(regular) singularity can be brought into the Fuchsian form. This means that if a G £ is 
a moderate singular point, then there exists a neighborhood U 3 a in CP 1 and a matrix 
function H(t) holomorphic and invertible in U\a having (together with P _1 (t)) at most 
a pole at a, such that the matrix 1-form dH ■ H^ 1 + HflH^ 1 has a Fuchsian singularity 
(first order pole) at a. This form is the logarithmic derivative dY ■ Y^ 1 of the matrix 
function Y(t) = H(t)X(t) ("linear change of the dependent variables"), provided that 
il = dX ■ X^ 1 is the logarithmic derivative of X. 

The problem of global classification of regular systems is closely related to the Hilbert 
21st Problem. It was discovered only relatively recently (by A. Bolibruch in 1990) that in 
general not every regular system can be brought to the Fuchsian form by a rational gauge 
transformation H with singularities only at the polar set: there is a subtle obstruction for 
systems with a reducible monodromy group. However, this is always possible if one allows 
creation of a single extra singular point outside £. Thus without loss of generality we 
may always assume that a regular system can be brought into the form @ by a rational 
gauge transformation. 

Remark 2.2. Together with the analytic description given above, the regular systems 
can be described geometrically. Given a regular system, with any simply connected open 
domain T a C CP 1 \ £ we can associate the n-dimensional linear subspace Jz? a in the 
space of analytic functions 0{T a ), spanned by solutions of the system. Any non-empty 
simply connected pairwise intersection T a @ defines the canonical isomorphism between 
J£ a and Jzfg, in other words, the spaces Jzf a are naturally globally organized into the holo- 
morphic vector bundle over CP 1 \ £ '. If the initial system is regular, then one can extend 
this bundle in an analytic way over the singular locus as explained in IY08 , Proposition 
18.8]. Moreover, the procedure of analytic continuation induces on this bundle the addi- 
tional structure of a connection in such a way that the functions from Jz? Q correspond to 
horizontal sections of this connection over T a . By construction, this connection is holo- 
morphic and flat over CP 1 \ £ . At the singular points the connection has an "atomic" 
curvature defined by the local monodromy operators. It is this geometric structure which 
naturally appears in many applications, which makes regular systems so ubiquitous in 
mathematics. 

2.2. Root counting for Fuchsian systems. The fact that the linear span of all entries 
of X is well defined, allows us to formulate the root counting problem for the system (JT]). 
In what follows we use the symbol #A for the number of isolated points in the set A. 
Consider the supremum 



sup sup 

T n £ = ca G c 




(3) 
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Here T C CP 1 \ S denotes an open triangle (an area bounded by 3 circular arcs or line 
segments) free from singularities of the differential equation (flj and supremum is taken 
over all such triangles free from singular points and all linear combinations ^ CyAy(i). 
If this supremum is finite, then we say that the system ((T|) admits an upper bound for the 
number of isolated zeros of solutions, or, simply, admits root counting. 

Example 2.3. A system with an irregular (wild) singular point t = a in general does not 
admit root counting. In fact, by general arguments from the theory of value distribution, 
any solution of super-polynomial growth admits almost every value infinitely often, in 
any sufficiently large sector around the singular point. 

For Fuchsian singularities the possibility of local root counting (say, in a small circular 
neighborhood of the singular point) depends on the spectrum of the residue of the matrix 
1-form O. 

Example 2.4 (continuation of Example 12. It is well known that if all the eigenvalues 
Xj of the residue matrix A — Aq = — A^ are real, then any such linear combination 
with real coefficients Cj may have no more than n — 1 positive real roots (the non- 
oscillation, or the Chebyshev property): one can prove this fact using the Rolle lemma 
on alternation between zeros of a real smooth function and its derivative. The complex 
roots of the linear combinations with complex coefficients can also be counted using a 
complex version of the Rolle lemma, as shown in [KY96J. In particular, under the same 
assumption Ai, . . . , \ n G M, the Euler system admits an upper bound for the number of 
roots by the expression (n — 1) + 2tt ■ |Ai — A n | (assuming the eigenvalues are labeled in 
a monotone order). 

On the other end, if n = 2 and Ai = — A2 = \f— 1 is a pair of complex conjugate 
eigenvalues, then the function ^(t 1 + i _i ) = cos In t has infinitely many positive isolated 
roots accumulating to both t = and t = +00. In this case N = +00 and the system 
does not admit root counting. 

In fact, the property that is essential for the root counting, is more naturally expressed 
in terms of the local monodromy, i.e., the monodromy operators associated with the 
small loops (around different singular points). For a Fuchsian singularity the eigenvalues 
Pi, . . . , Vj of the monodromy of a small loop are exponents, Vj — exp 2iri\j, of eigenvalues 
Ai, . . . , A„ of the residue matrix |IY08[ Corollary 16.20]. Thus all eigenvalues of the residue 
are real if and only if all eigenvalues of the local monodromy operators have unit modulus, 
\vj\ = 1. (Note that this description fails for regular non- Fuchsian singularities). 

Definition 2.5. We say that a regular system (fTJ) has quasiunipotent monodromy (or 
simply is quasiunipotent), if the monodromy operators Mj of all small loops around 
distinct singular points aj € S are roots of unity (in particular, have only modulus 1 
eigenvalues). 

Using the technique developed by A. Gabrielov and A. Khovanskii, one can prove 
the following (relatively easy) theorem, whose assumptions are algebraic and can be 
effectively verified. 
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Theorem 2.6. A regular Pfaffian system ([T} with quasiunipotent local monodromy admits 
finite root count, N(f2) < +00. 

Unfortunately, this approach, albeit very general and powerful, does not allow to make 
the bound for the counting function explicit, even for Fuchsian systems. 

2.3. Explicit bounds for Fuchsian systems. Several simple examples show that al- 
though each quasiunipotent Fuchsian system ([2]) with the Pfaffian matrix SI having only 
simple poles with real spectrum residues, admits a finite bound N(f2), this bound cannot 
be uniform on all systems of the given degree m (the number of singular points) and 
dimension n. 

Example 2.7. The diagonal Euler (2 x 2)-system with two integer eigenvalues and 
I 6 N admits the complex polynomial t — 1 as a solution; this polynomial has as many 
as £ isolated roots in a suitable triangle T C CP 1 . This suggests that the bound N(f2) 
may well grow to infinity as the spectral radii of the residues grow to infinity. 

A slightly less obvious obstruction to uniformity of the bound is the sensitivity of the 
quasiunipotence condition to the "collision" of singularities, when the distance |<Zj — a,j\ 
between two distinct singular points ai, aj g S tends to zero. It may well happen [BY09 
that the monodromy along each small loop is quasiunipotent, but the loop which encircles 
both singular points, has non-quasiunipotent monodromy. Thus in the limit when aj = aj, 
the system may have infinitely many roots of solutions, and an arbitrarily large finite 
number of them would persist for < \di — a,j\ <C 1. 

These phenomena can be re-stated as follows. The collection m of all Fuchsian sys- 
tems ^) of dimension n with m singular points is a semialgebraic variety, a subset of a 
complex affine space which can be represented as the difference of two affine algebraic 
varieties. For instance, consider all possible residue matrices A%, ... , A m e Mat(n, C) ~ 
C" and all possible pole assignments ax,..., a m <= C, and exclude the "degenerate cases" 
when Ai = (disappearance of one of the singular points) or aj = aj (collision of the 
poles). Then the counting function N(-) becomes well defined on a proper semialgebraic 
subset § n ,m 5 3>s,m defined by the condition Spec(A,) CR,j = l,...,m. The counting 
function 3Sf is lower semicontinuous (all isolated complex roots of linear combinations 
survive sufficiently small perturbations if counted in open triangles with proper multi- 
plicities) and hence this function is bounded on any compact subset of §> n ,m- However, 
since §„. m itself is non-compact, N(-) can apriori be unbounded, and the above exam- 
ples show that indeed N tends to infinity at least near some parts of the boundarj0 

The main result of the paper [BY09] is the claim that the counting function N grows no 
faster than polynomially near the boundary dS„ tm . This means that for any semialgebraic 
"reciprocal distance to the boundary" p(fl) = l/dist(f2, dS ntm ) one has an estimate of 
the form 3ST(fi) ^ p k (ft), where k < +00 is a finite power which depends, of course, on 
the choice of p (and also implicitly on n and m). Since by the Lojasiewicz inequality any 

1 More precisely, we have to consider the projective boundary for which the closure should be 
taken with respect to the projective space which compactifies the affine space C . 
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two such distances p, p' are equivalent, p' p°^ and vice versa, the above polynomial 
growth statement is in fact independent of the choice of the distance p as long as the 
latter remains semialgebraic. 

However, the real strength of the bound achieved in |BY09j is its constructive nature: 
it allows to give an explicit upper bound for k in terms of n,m, assuming any given 
distance function. For instance, for a Pfaffian system ([T]) with residue matrices Aj and 
the singularities aj (in a given affine chart) we can define 

p (n) = 2 + J2\A j \ + J2r^-\ <+(X3 - W 

(In fact, this function does not "notice" the disappearance of the singular points when 
Aj —> 0; the constant 2 is added in order to avoid special treatment of systems with 
p ^ 1, as we are interested only in large values of p). 

Theorem 2.8 ( [BY09] ). 

nn)^ P {n) 2 \ (5) 

where Poly(n, m) is an explicit polynomial expression in n, m of degree at most 1 + 
dimS„, m = l + mn 2 + m, where dimS„. m is the dimension of the semialgebraic setS niTn . 

This double exponential bound may seem too excessive, yet one can argue following 
[BY09] that a non-uniform bound polynomially growing near the boundary of a semial- 
gebraic set in a high-dimensional affine space and invariant by the choice of the semial- 
gebraic distance, should necessarily be growing at a comparable rate as the dimension of 
the space grows to infinity. 

2.4. Towards uniform bounds: obstructions and solutions. As was already men- 
tioned in the previous section, one cannot expect to obtain a bound which would be 
uniform over all domain of finiteness of the counting function N. Thus one has to focus 
on constructing some smaller "naturally noncompact" subsets of S„. m , for which uniform 
bounds can be achieved. 

The weakest (also the most important) barrier occurs on the diagonal set Ui^j{ a i — 
aj}. As was explained above, the obstruction to existence of a uniform bound is hidden in 
the fact that after collision of two singular points with quasiunipotent local monodromy 
one can obtain a point with non-quasiunipotent monodromy, to which infinitely many 
roots may accumulate. To prevent this, one could impose the quasiunipotence condition 
also on the monodromy for the loop encompassing both colliding singularities. However, 
this condition is non-algebraic and in general non- verifiable: unlike the local computation 
with finite order jets which allows to compute the local monodromy operator [IY08, 
Corollary 16.19], the monodromy along a large loop requires in general knowing exact 
solutions which are generically non-algebraic. 

Instead one can impose an additional condition that all systems from the conjectural 
class are isomonodromic to each other. This condition stipulates that for each closed 
loop avoiding the singular locus of a given system, all sufficiently close systems from this 
class have the same monodromy along this loop modulo conjugacy. In particular, this 
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condition implies that all residue matrices remain in the same conjugacy class and hence 
their characteristic polynomials are constant. 

The isomonodromy condition is surprisingly rigid: in general, it (locally) uniquely 
determines the residue matrices as functions of the polar configuration. More precisely, 
in the so called non-resonant cased the residues A4 — Aj(ai, . . . , a n ) of the Fuchsian 
system ([2]) considered as functions of the pole location 01, . . . , a n , satisfy the integrable 
system of (nonlinear) matrix Pfaffian equations, called the Schlesinger system: 

dA t + VL4„ Aj] d( " Ql ~ a ^ = 0, i = l,...,n. 

a,i — aj 

Solutions of this system may exhibit singularities both on the diagonal ( "collision locus" 
Ui=y{ a « — a j — 0}) an d elsewhere (depending on the initial conditions), since the system 
is nonlinear (quadratic). Therefore in general the bound ^ will diverge because of the 
"explosion" of the residues in the isomonodromic collision of Fuchsian singularities. 

This phenomenon suggests that parametrization of isomonodromic families of Fuch- 
sian systems by the location of their poles and the corresponding residues leads to creation 
of singularities on the diagonal (the collision locus) . The alternative is to change the point 
of view and consider parametric classes of regular linear systems as (Pfaffian) systems on 
multidimensional phase space. Then the Schlesinger equations can be interpreted as the 
holonomy (flatness) condition equivalent to the local existence of solutions. This point of 
view is formalized as follows. 

Consider a projective space CP m and a rational (n x n)-matrix Pfaffian form fl = 
{f2jj}™ J=1 on it. Denote by 2J C CP m the polar locus of fi, the union of polar loci of all 
scalar forms fi^ € /\ 1 s (CP m ). Outside of 2J we can consider the system formally defined 
by the same equations as ([T|) , yet this time with several "independent variables" , 

dX = nX, ft e Mat nx „(A^(CP m )), X € GL(n, ff{CP m \ U)), (6) 

which becomes a system of partial (rather than ordinary) differential equations in an 
affine chart on CP™. Obviously, we are interested only in systems which locally admit a 
holomorphic fundamental solution matrix X so that f2 = dX ■ X^ 1 . Applying the exterior 
derivative to both parts of this formula, we conclude that 

dfi = -dXAdiX- 1 ) = -dXAi-X^dX-X- 1 ) = (dX-X" 1 )A(dX-X^ 1 ) = flAfl. (7) 

This local integrability condition is also sufficient for the local existence of solutions; it 
has a geometric meaning that the connection defined by the form Q on the trivial bundle 
over CP™ is flat, i.e., has zero curvature. 

In the same way as in one dimension, we define the regularity condition for the system 
(j6]) on the singular locus, by requiring that any solution grows no faster than polynomially 
when approaching £ inside any simply connected semialgebraic set U C CP™ \ S (the 
semialgebraicity condition is required to exclude sets spiraling around S): 

\X{z)\ + < Cdist{z,S)- N for some finite N < +00, z G U. (8) 

2 The non-resonance assumption requires that the difference between any two eigenvalues of 
each residue matrix Aj is non-integer. 
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The root counting problem for solutions of the system © (defined on CP 1 ™) by defini- 
tion reduces to that for the system ([l) (defined on CP 1 ) by considering the restriction on 
all possible projective lines not lying entirely in the singular locus: formally we can recycle 
the definition ([3]) by adding that the supremum is taken over all triangles T belonging to 
all projective lines in CP 1 . 

The next section contains the sufficient conditions for a regular integrable system |6]) 
on CP m to admit finite root count. We formulate a proper multidimensional generaliza- 
tion for the quasiunipotence condition on the local monodromy operators of a regular 
system ((T|), which guarantees (together with the above conditions of integrability and 
regularity) the finiteness of roots in any triangle. 

2.5. The main result on quasiunipotent systems. 

Definition 2.9. A small loop around the polar divisor £, centered at a point a £ £, 
is the free homotopy class of the image of any sufficiently small circle {\t\ — e} by a 
holomorphic map 7 : (C, 0) — > (CP m , £) whose image does not belong to £. 

If a € £ is a smooth point of £, i.e., £ is locally near a represented by a single 
equation {/ = 0} with f(a) = 0, d/(a) 7^ 0, then all small loops centered at a are 
free homotopy equivalent to each other. On the other hand, if a is a singular point of 
the locus £, this may not be the case anymore. However, certain spectral properties of 
the corresponding monodromy operators are inherited from generic to degenerate small 
loops. 

Theorem 2.10 (Kashiwara |Ka81j ). Assume a flat regular system © with the polar locus 
£ possesses the following property: its monodromy operator along any small loop centered 
at a smooth point £ is quasiunipotent. 

Then the same quasiunipotence hold for any small loop, including those centered at 
the singular points of £ as well. 

This theorem is in fact a result on the local topology of the complements to embedded 
analytic hypersurfaces. Note that all small loops centered at smooth points of £ are also 
free homotopic to each other in the case when £ is an irreducible algebraic subvariety 
in CP™. Thus the quasiunipotence condition is sufficient to verify at a finite number of 
points (one smooth point for each irreducible component of £). Moreover, if f2 has a first 
order pole on £ (i.e., / f2 is holomorphically extendable for any local equation of £), then 
one can correctly define the residue as a holomorphic matrix function on the smooth part 
of £, such that all values of the residue are conjugate to each other along any irreducible 
component of £. For such systems the quasiunipotence condition is equivalent to the 
assumption that the residue matrix has only rational eigenvalues. 

A system © meeting the assumptions of the Kashiwara theorem, will be called qua- 
siunipotent, naturally extending thus Definition 12.51 Not surprisingly, after the notions 
of regularity ((H)), quasiunipotence and the root counting function N are generalized from 
the one-dimensional to the multidimensional case, Theorem 12.61 remains valid. 

Theorem 2.11. A quasiunipotent flat regular Pfafflan system © admits finite root 
count: N(f2) < +00. 
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As its one-dimensional version, the mere finiteness of the bound can be proved using 
very general arguments of Gabrielov-Khovanskii-Varchenko type. To make this bound 
explicit, a completely new approach is required. Even to formulate such a result, one has to 
introduce an additional characteristic of linear systems (O which would essentially restrict 
the "magnitude of the coefficients" to exclude the situation described in Example 12.71 
Such a characteristic must necessarily enter the expression of any explicit bound, and the 
discussion in £12.31 suggests that the norm of the residue matrices is a wrong choice. 

Definition 2.12. We say that a Pfaffian system is defined over Q, if in some affinc 
chart C' m = {ti, . . . , t m } on CP™ all rational 1-forms ttij have coefficients from the field 
Q(ti, • • • , t m ) 1 i.e., ratios of polynomials with integer coefficients from Z[ti, . . . , t m ]: 

m P i (t) 

The largest integer number which is used to write explicitly the polynomials Pijk,Qijk, 
is called the complexity of a system defined over Q. 

This definition clearly depends on the choice of the representation of (e.g., it changes 
after the affine change of variables 1 1— > 2t) , yet for our purposes it will always be used in 
the sense that at least one representation of the given complexity exists. Having an upper 
bound on the complexity implies that any numeric characteristic of the system, that can 
be obtained by a finite number of algebraic manipulations, can be explicitly majorized in 
the sense of the absolute value. This allows to translate any explicit algebraic algorithm 
into an explicit bound for the complexity of its results. In particular, we obtain the 
following bound for the root counting problem. 

Theorem 2.13 ([ BNYIO I). If a quasiunipotent flat regular Pfaffian system © is defined 
over Q and has complexity at most s £ N, s ) 2, then 

W(fi) < s 2 Pol y( rf ' n ' m ) ) Poly(d, n,m) < 0((d 5 m 5 n 20 )). 

Here d = maxjj^deg Pijk, deg Qijk) is the degree of the Pfaffian form, m — dimcCP m 
the dimension of the space [number of the independent variables) and n — dimfi the 
dimension of the system (number of the dependent variables). The polynomial of degree 
30 occurring in the exponent is explicit and can be computed. 

This theorem implicitly introduces a class of special "quasialgebraic" multivalued 
functions, which generalizes the class of algebraic functions in the sense that they admit 
finite root count of any branch. 

Definition 2.14. A Q-system is a flat regular quasiunipotent Pfaffian system defined 
over Q. A Q-function is a multivalued analytic function which is a component (or linear 
combination of components) of a matrix solution of a Q-system. 

The assertion that a multivalued nondegenerate matrix function X(t), t E CP m \ S, 
is a Q-function, consists therefore of four separate statements: 

1. X undergoes a linear monodromy transformation X i-> XM, M 6 GL(n,C), when 
continued along a loop 7 avoiding the ramification locus S; 
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2. This monodromy transformation has only roots of unity as the eigenvalues, if 7 is 
a small loop (by the Kashiwara theorem, this condition needs to be checked only 
for small loops around smooth part of S, one loop for each irreducible component 
of S suffices); 

3. X is regular, i.e., exhibits at most moderate growth on £ (in particular, X has 
genuine finite order poles along the components of S with the identical monodromy 
or satisfies the growth control estimate ([5])): 

4. The logarithmic derivative — dX ■ X , necessarily rational in t under the above 
assumptions, is defined over Q. 

The flatness condition is automatically satisfied. 

3. Q-systems, Q-functions and their properties. 

3.1. Elementary examples of Q-functions. Polynomials C[tx, . . . , t m ] of a given 
bounded degree d are a trivial example of Q-functions. Indeed, the exterior derivative 
of each monomial t v = t v ^ ■ ■ ■ is a linear combination of monomial 1-forms of smaller 
degree with natural coefficients not exceeding d, which means that the corresponding 
logarithmic derivative in this affine chart is constant and defined over Z. The growth 
condition is obviously satisfied and the monodromy is identical, hence trivially quasi- 
unipotent. Of course, the double exponential root count provided by Theorem 12.131 is 
enormously excessive in this case. 

Proposition 3.1. Any algebraic function defined over Q is a Q-function. 

Before proving the Proposition, we need a simple algebraic fact which is a cornerstone 
of the Elimination theory. 

Lemma 3.2. Let S c C" x C be an algebraic hypersurface defined by the polynomial 
equation S(t,x) = of degree d — deg x S, with S defined over Q. Then any rational 
function R(t,x) = A(t,x)/B(t,x) € Q_(t, x) restricted on S coincides with a fraction of 
the special form, 

U(t x) 

R(t,x)= frj ' U e 1[t,x), Q e Z[t], deg x U^d-l. 

The new denominator Q depends only on the polynomial S defining the surface and the 
denominator B of the rational function R. The complexity of this representation can also 
be explicitly controlled in terms of the complexities of R, S and their degrees. 

Proof of the Lemma. Consider the resultant of S (the equation) and B (the denominator) 
considered as polynomials in x with coefficients in the field Q(t). By definition, this is 
the element of the field Q(t) representable as the linear combination, 

Q(t) = a(t, x)S(t, x) + b(t, x)B(t, x), a,b e Q(t, x), 

with deg^, a < deg x B, deg x b < deg x S. The resultant Q(t) vanishes if and only if S(t, ■) 
and B{t, •) have a common x-root. 

Evaluating this identity on the hypersurface 5*, we conclude that 1/B equals the poly- 
nomial b/Q. Dividing A(t, x)b{t, x) with remainder by S(t, x) (using the Euclid algorithm) 
in the ring Q(x)[t], we find U, V e Q(x)[t] such that Ab = VS + U with the required 
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bound for the degree deg^ U ^ deg^, S. Multiplying the numerator and denominator by 
a suitable term from Z[x], we can guarantee that the ratio involves polynomials in t,x 
with integer coefficients. ■ 

Proof of the Proposition. The algebraic equation P{t\, . . . , t n , y) = defines an algebraic 
surface in C™ x C, and its exterior derivative dP = implies the Pfaffian equation 

d y = -i^7r^J2 Pkdt ^ p, = ^r> Pfe = ^ p',p k eQ[t, y \. 

Applying Lemma I3.2[ we can replace the restriction of this equation on P = by a 
rational Pfaffian equation whose denominator A(a;), the resultant of P(t,y) and P'(t,y) 
as polynomials in y (the discriminant of P(t, •)), depends only on t and the numerator 
has degree in y not exceeding d — 1 : 

1 ™ 

^^E^t^)*' U k €Z[t,y], deg y U k ^d-l. (9) 

^ ' k=l 

The same procedure can be applied to the Pfaffian equation @ multiplied by any power 
y 1 . the result will be a rational Pfaffian form for the exterior derivative d(y^), j = 
2,3, .. . whose denominator is still the same discriminant A(t) and the numerators are 
polynomials of degree ^ d — 1 in y. 

In other words, we obtain a linear system of Pfaffian equations 

_^ d n 

dXj = Mi ^ ^ u jkl (t)x k dti 
' ' fc=i i=i 

for the variables x%, . . . ,x n , xj = y^~ x , j — 1,. . . ,n, with coefficients from Q(t) (obvi- 
ously, dxi = 0). The algorithm is explicit and allows for estimation of the corresponding 
complexity. 

The monodromy of this linear systems simply permutes between themselves the 
branches of the algebraic function. Any such transformation is cyclic (its Jordan blocks 
are cycles of period no greater than n) hence the corresponding monodromy is quasiu- 
nipotent for any loop. ■ 

Corollary 3.3. The general algebraic function defined by a "universal polynomial equa- 
tion of degree n " , 

to y n + h y n - x + ■■■+ t n -i y + t n = 0, (10) 

is a Q-function of its arguments to,---,t n considered as homogeneous coordinates on the 
projective space CP™. ■ 

3.2. Special functions. Hypergeometric and Riemann P-functions also are Q-functions 
for the rational values of the parameters almost by definition: the second order linear 
equation 

d^?y dv 
t(t — 1) • -777 + [c — (a + 6 + l)t] -£-ab-y = 
dt z dt 

has three regular singular points at t = 0, 1, 00 and the corresponding residues are rational 

for rational a, b, c € Q. Thus a vast number of special functions turns out to be Q-functions 

and admit explicit bounds on the total number of roots (depending on a, b, c). 
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However, as a function of all four variables t, a, b, c, the Ricmann P-function is not 
a Q-function: the parametric family above is not isomonodromic and hence cannot be 
transformed to an integrable Pfaffian system on CP 4 . Besides, some solutions grow faster 
than polynomially as (a, b, c) — > oo. 

3.3. Periods as Q-functions. The first general example of transcendental Q-functions 
is provided by Abelian integrals 

Definition 3.4. An Abelian integral is the period /(P, w) of a rational 1-form oj e 
/\ (CP 2 ) over an algebraic cycle on an algebraic curve P C CP 2 , considered as the 
function of all relevant data. 

Recall that an algebraic cycle is an element c of the homology group Hi (P, Z) of 
an algebraic curve P C CP 2 . If the curve is nonsingular, then this cycle (as a locally 
constant section of the Gauss-Manin connection) continuously depends on the curve P 
in a natural way. In an affinc chart (x, y) on CP 2 an Abelian integral takes the form 



The function / is naturally defined on an the Zariski open subset of nonsingular algebraic 
curves of a given degree n, parameterized by points of the suitable projective space 
(coefficients of the polynomial H) and the projective space of rational forms of a given 
degree ^ m. 

Theorem 3.5. The Abelian integral (|llj) is a Q-function of its arguments. The dimen- 
sions, degree and complexity of this Q-function are bounded by a polynomial in the degrees 



Remark 3.6. The fact that the Abelian integrals are regular functions satisfying some 
linear (Picard-Fuchs) ordinary differential equations with rational coefficients, is well 
known under much more general assumptions. Consider an arbitrary surjective polyno- 
mial map P : C" — > C m , n ^ m, denote by T t = P _1 (£) the fibers of this map, which are 
affine algebraic varieties. For a generic value of t (i.e., off an algebraic surface £ C C m ) 
Pt is nonsingular variety, whose homology is generated by algebraic cycles of dimension 
k = n — m which locally continuously depend on t ^ S. By a theorem of Grothendieck 
|Gr66], the dual module (cohomology) is generated by the same number of polynomial 
fc-forms whose restrictions on P t are linear independent outside a proper polynomial hy- 
persurface (we can still denote it by £). Thus the period matrix X(t), whose entries are 
all periods of the corresponding forms, is well-defined, holomorphic and nondegenerate 
outside £ matrix function, which grows moderately as t — > £ and is monodromic: after 
analytic continuation along any loop 7 around S, the matrix function acquires a right 
(constant) matrix factor M 1 completely determined by the topology of P considered as 
a bundle. The quasiunipotence of M 1 is also a well-known fact. The matrix logarithmic 
derivative ft = dX ■ X^ 1 is a rational matrix 1-form on C m (hence extends on CP m ) 
with poles on £, and the first three conditions from Definition 12 . 141 are verified. 

However, this quite general argument does not allow to see why the matrix f2 is defined 
over Q, and does not allow to estimate its complexity. One has to construct explicitly the 





ofP,Q,Z,H. 
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matrix fi in such a way that its complexity over Q can be bounded. 

This construction was produced in [BNY10 for the particular case where the polar 
locus {Z = 0} of the forms is a fixed line in CP 2 so that the space of all forms is a 
linear space of polynomial 1-forms in a fixed affine chart. The general case will be treated 
elsewhere. 

3.4. Elementary constructions with Q-functions. The ring operations (addition 
and multiplication) preserve the class of Q-functions, affecting only their natural charac- 
teristics (degree, dimension, number of variables and complexity). 

Indeed, for two Q-systems defined on two different projective spaces by two systems 
of quasiunipotent flat regular Pfaffian equations 

dx = nx, dY = ey, 

of the sizes n x n and m x m respectively, one can form the block diagonal system of size 
(n + to) x (n + to), 



d 



which is obviously flat, regular and quasiunipotent (the monodromy group is the Carte- 
sian product of the two individual groups). This system (whose matrix can be denoted 
by fl © 0) is satisfied by the direct sum X © Y of the corresponding matrix solutions. 

The tensor product X®Y, also satisfies the linear system of the size nm x nm. To see 
this, it is more convenient to look at the two systems in the vector (rather than matrix) 
forms: 

dx = 0.x, dy = VLy, x = (xt,. .. ,x n ), y = (j/i, . . . , y m ). 

Then the n x m-vector with the coordinates = x^yj satisfies the system of linear 
Pfaffian equations 

n m 

dzij = dx t ■ y 3 + x t ■ dy } = ^ ClikZkj + QjiZn- 

k=l 1=1 

The matrix elements of the corresponding rtm-matrix 1-form are certain sums of the 
matrix elements of the forms £1 and respectively. We abbreviate the corresponding 
computation using the tensor product notation as follows^), 

d(X ®Y) = dX®Y + X®dY = (rL®I + I® G)(X <g> Y). 

One can easily verify that all characteristic properties of the Q-systems are preserved by 
the tensor product, with explicit control over the complexity growth. 

These elementary properties imply that the extension k(X) = C(i, X) = C(ti, . . . , t m , 
Xu, . . . ,X nn ) of the differential field k = C(ti, . . . ,t m ) by Q-functions X = X(t) is a 



3 In the language of connections one sometimes uses the following notation for the tensor 
product, 

(V' <g> V)(X <g> Y) = (V'X) ® Y + X ® (V"F), V' = d - n, V" = d - 6. 

This may justify the notation Q ® Q for the matrix form corresponding to the tensor product 
V <g> V" — d — Q (g) O of the two connections. 
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differential field which admits an explicit upper bound for the number of isolated roots 
of each its element. 

Theorem 3.7. Let X be any fundamental matrix solution of a fixed Q-system © and 
k.(X) the differential extension field. 

Then for any element <p € h(X) and any one- dimensional triangle T entirely belonging 
to the domain of analyticity of if, the number of isolated zeros of if in T is bounded by 
an explicit function of the degree deg t x f '■ 

VP c CP m \S, V^E k(X) ff{t e T : <p(t) = 0} < 2 2P ° IyW , S = deg tjX <p. (12) 

The polynomial Poly((5) explicitly depends on the Q-function X via the dimension, degree 
and complexity of the latter. 

Proof. All monomials of the form t a X^ with \a\ + \j3\ ^ 5 (in the standard multi-index 
notation) satisfy a large Q-system whose dimension, degree and complexity are bounded 
by some explicit polynomial expressions in the dimension, degree and complexity of the 
Q-function X, and the degree S. The rest follows from the Main Theorem 12.131 ■ 

3.5. Transformations of Q-functions by rational and algebraic maps. If X is a 

matrix Q-function on CP m and / : CP' — > CP m a rational map defined over Q, then the 
natural pull-back f*X will be a Q-function on CP', defined by the Pfaffian system with 
the matrix 1-form /*f2. 

Indeed, the conditions of flatness and regularity are obviously satisfied. The pullback 
matrix 1-form f*Q is clearly defined over Q. The only condition that has to be verified 
is that on the monodromy. But the image of a small loop by a rational map is at worst 
a multiple of a small loop, hence the eigenvalues of the corresponding monodromy are 
powers of the initial eigenvalues. 

Yet a much more interesting question is to find out under what conditions the class 
of Q-functions is closed by composition, in particular, when the substitution of the in- 
dependent variables t = F(s) with algebraic functions F transforms Q-functions again 
to Q-functions. Among other things, this would allow to replace "line triangles" in the 
formulation of the counting problem, by arbitrary semialgebraic triangles. 

Theorem 3.8. Consider a Q-system (O and a transformation defined by an algebraic 
change of independent variables, 

tj =tj(si,...,si), j = 1, P j (t j ,s 1 ,...,si) = 0, PjEQ[tj,s]. (13) 

Then the multivalued matrix function X(t(s)) is a Q-function of the new variables. 

Proof. Consider the Pfaffian system ([6|) in the affine chart (ti, . . . , t m ): 

dx t = 22Rijk{t)xjdt k , Rij k €Q(h,...,t m ), degPyfc^d. 

Differentiating the identities (fT3"|) . the differentials dtj can be replaced by rational forms, 

ft 

After substitution we obtain the linear Pfaffian system with the "dependent variables" 
xi, . . . ,x n and independent variables s%, . . . , s;, although the coefficients are not rational 
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anymore: they explicitly involve the algebraic functions tj(s). However, this system can 
be transformed to a system with coefficients in Q(s) in a way similar to that used in the 
proof of Proposition 13. f I 

Denote by k = Q(s) the field of rational functions of the variables s. Reducing all 
equations to the common denominator, we can transform the Pfaffian system to the form 



with the coefficients Qijk polynomial in t. Without loss of generality, we may assume that 
all such polynomials are of degree ^ d — 1 by virtue of the identities (|T3^) . 

We first replace the denominator A(t) by an element from the field k = Q(s), by iter- 
ated application of Lemma 13.21 More specifically, we eliminate inductively the variables 
t±, . . . , ifc from the denominator, replacing them by the variables Sj, using the equations 
(|13|) . Then the standard division with remainder allows to exclude all powers of tj of 
degrees higher than dj — deg tj Pj. Finally, the lower degree powers tj,tp . . . 1 are 
declared to be the "new dependent variables" governed by the obvious recurrent equations 
d(fj) = rt T j~ 1 dtj transformed to the linear Pfaffian form using the identities (HU). 

The quasiunipotence the monodromy is proved as follows. Consider a small loop 
7 : (C,0) — > (CP', a), defined by a holomorphic germ s — s(z), and its composition 
with the algebraic map r\ : z > t = t{s{z)). Since the algebraic maps are in general 
multivalued, the latter map in general is multivalued (ramified at the origin) and its 
restriction on a small circle \z\ = p non-closed. However, the algebraic map (fT3"|) has only 
finitely many branches which are permuted along the small loop 7. Thus after finitely 
many iterations the small loop 7^ would become a closed loop in the t-space. On the 
level of the monodromy operators this means that the monodromy Af 7 of the composite 
system along the small loop 7 is quasiunipotent when raised to some finite power TV, i.e., 
all eigenvalues of are roots of unity. But then the same is true for the eigenvalues of 
M 1 itself. ■ 

Remark 3.9. Inspection of the above proof shows that one can use a more general 
transformation of a "triangular" form, where the coordinate functions of the algebraic 
transformation are defined by the "Pfaffian chain" of the polynomial equations 



P(tx,s) = 0, P 2 (h,t 2 ,s)=0, ... P m (t 1 ,...,t m ,s)=0, 8 = («!,..., sj)- (15) 



It would be very interesting to know whether the restriction of a Q-system on an arbitrary 
projective subvariety Z defined over Q, can always be represented by a Q-system (with 
eventually larger number of the dependent variables). On the other hand, probably there 
are compositions of Q-functions, which are themselves not Q-functions. 

Remark 3.10. The possibility of introducing extra parameters and additional auxiliary 
equations (algebraic or Pfaffian) implies that in the definition of a Q-system one can 
replace the assumption that Q is defined over Q by formally weaker but in fact equivalent 
condition that it is defined over the ring of algebraic numbers Q. 





4. Applications to the Hilbert 16th problem. Although the above exposition was 
aimed at showing that the class of Q-functions is rich enough and important by itself to 
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merit an investigation, still by far the strongest motivation for its study comes from the 
theory of planar real algebraic foliations and bifurcations of limit cycles. We recall briefly 
the connection between these two areas. Detailed exposition can be found in numerous 
textbooks, among them |IY08j . 

4.1. Crash course on perturbations of Hamiltonian foliations. If H € M[a;,?/] 
is a real polynomial in two variables with isolated critical points, the Pfaffian equation 
dH = defines a real foliation with singularities on the plane M 2 which extends to the 
projective plane M.P 2 . Leaves of this foliation are connected components of the real level 
curves {H = t}. This foliation is integrable: the holonomy associated with each compact 
leaf (real algebraic oval) is trivial (identical), in particular, all nearby leaves are also 
compact and closed. In the classical language of the differential equations, there are no 
limit cycles for the corresponding system of Hamiltonian equations. 

This changes if we apply a small one-parameter perturbation and consider the Pfaffian 
equation 



dH + euj = 0, uj = Pdx + Qdy, P, Q, H G R[x, y], £€^,0). (16) 



Since the holonomy of a compact leaf is analytically dependent on the parameter e, it 
can be expanded as a converging series. Using the function H as the local chart on the 
transversal, we can write 



where Ii(t), l2(t), ■ ■ ■ are the first and the higher variations of the holonomy, functions, 
defined on the nonsingular compact leaves of the initial integrable foliation and analyti- 
cally depending on these leaves. A simple calculation known already to Poincare shows 
that 



is the real Abelian integral, cf. with (jTTJ) . If the first variation is nontrivial, h(t) ^ 0, 
then the compact isolated leaves (limit cycles) of the perturbed foliation (|TH|) are rather 
faithfully tracked by the real isolated roots of the period (fTB|) (an accurate statement 
involves the positive distance from the singular locus of the integrable foliation). If I\ = 0, 
then one should find the lowest order nontrivial variation k ^ 2, and study its zeros. 

This construction justifies the following "linearized" ("infinitesimal") relaxed version 
of the Hilbert's 16th problem: find an upper bound for the number of isolated zeros of 
the Abelian integral (TT51) in terms of the degrees degH and degw = max(deg P, deg Q). 

Clearly, the role of the two ingredients, the Hamiltonian and the perturbation form, is 
different: e.g., as a function of the form, the integral (|18j) is linear, while its dependence on 
H is very tricky. As an intermediate step, one can split the problem, fix H and consider 
the problem for various 1-forms uj of different degrees. 

One can show that for each fixed H of degree n one can find finitely many polynomial 
1-forms which generate all integrals (|18p as a module over the ring C[<] of polynomials: 
for a generic H one can take (n — l) 2 monomial forms u>i whose differentials du>i span 
the finite-dimensional quotient space /\ (R 2 )/ A (^ 2 ) A di?„, where H n is the principal 
homogeneous part of H. The assertion means that for any polynomial 1-form u of degree 



rj(t,e) = t + sh{t)+s 2 I 2 {t) + -- - , 



(17) 





(18) 
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d and any continuous family of cycles Ct C {H = t} there exist polynomials qj G M[t] 
such that 

j> oj = Y,<lj(t) ■ f e R[i], deg 9i sC d/n. (19) 

This means that the Abelian integrals of all polynomial 1- forms belong to the Picard- 
Vessiot differential extension field t(X), where X is the period matrix of the (n— l) 2 -tuple 
of 1-forms ujj over all cycles on the complex algebraic curves {H = t} C C 2 . 

4.2. Linear bounds. Theorem 13.51 allows to produce an upper bound for the number 
of isolated roots of a combination (|19p which would be double exponential in d. However, 
for quite some time it was known that this number admits an asymptotic estimate which 
is linear in d: it does not exceed C(n)d + C'(n), see [P97] and |Z061 Theorem 6.26, 
pp. 177-183]. Here the "constant" C(n), which depends only on the degree n — degH, is 
completely explicit, but the other constant C'(n) is only proved to exist, with the proof 
giving absolutely no clue on how it can be estimated. 

It turns out that this is a general feature of Q-functions. Consider a Q-system ((SJ on 
CP m and the module ~#(X) its solution matrix X{t) generates over the field C(f i, . . . , t m ) 
of rational functions on the domain of the Q-system: by definition, it consists of all 
linear combinations • tfoj (t)Xij (t) with rational coefficients qij(t) — q^ (t\, . . . , t m ). 
This representation, albeit not unique, induces a filtration of this module by the degrees 
of the coefficients g 3 : we will say that an element / € ^0(X) has degree 6 = deg t /, if it 
can be represented as a linear combination as above with deg gtj ^ 5 for all i,j. 

Theorem 4.1 f |BDll| . to appear). For any fixed Q-system of degree d dimension n and 
complexity s on CP m , the number of isolated zeros of any element of degree S from the 
module ^V(X) is bounded by an explicit expression 

C(n,m,d)-5 + C'(n,m,d,s). (20) 
Here C, C are two explicit functions bounded as follows: the leading coefficient C(n, m, d) 
grows no faster than a double exponential 2 independent of n, while the constant 

term C'(n,m,d, s) is bounded by s N ( n > m > d ) with the degree N bounded by a tower of five 
exponentials in (n,m,d). 

The proof is achieved by combination of the construction of "folding" suggested by 
G. Petrov and further elaborated by A. Khovanskii, with the effective bound provided 
by Theorem l3. 81 above. It gives a fully constructive linear bound for the number of limit 
cycles, born by polynomial perturbations of any degree, from a given polynomial Hamil- 
tonian foliation, as explained in Needless to say, this bound is also tremendously 
excessive by all expectations. 

4.3. Iterated integrals. The higher variations Ik(t) in (|T7|) can be computed recur- 
sively, see }Y95| IFr96j . For a generic Hamiltonian H it can be shown that, under the 
inductive assumption that I\(t) = ■ ■ ■ = Ik-i(t) = 0, the function /&(£) can be expressed 
as an Abelian integral of a polynomial 1-form of degree growing with k |G98| . see also 
}IY08j . Thus the number of limit cycles that can be tracked using variations of any given 
order k (and, as before, distant from the singular locus) , can be estimated in terms of k 
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by virtue of Theorem 14. II However, the problem of determining the maximal order k such 
that the identical vanishing of all I± , . . . , Ik implies that the perturbation (|16[) entirely 
consists of integrable systems, is a transcendentally difficult problem which includes as a 
particular case the famous Poincare problem of discrimination between center and focus. 

In the degenerate cases the higher variations Ik may not be periods of polynomials 
forms ( Abelian integrals) . For instance, the monodromy operators may have Jordan cells 
of size greater than 2 (which is impossible for the "genuine" integrals over cycles of 
dimension 1), see |GI05j . Nevertheless, one can always express the first nonvanishing 
variation Ik using iterated path integrals |G05j . These are special functions which in some 
sense interpolate between (dual spaces to) finite-dimensional homology of the complex 
fibers -H" -1 (£) and their infinite-dimensional homotopy. L. Gavrilov in [G05 proved that 
the iterated integrals considered as functions of a single variable t satisfy a regular system 
of linear ordinary differential equations. The monodromy of this system was shown to be 
quasiunipotent in [GNlOj . The explicit derivation of this system, achieved by S. Benditkis 
and D. Novikov, shows that this system is defined over Q. 

Theorem 4.2 (S. Benditkis and D. Novikov [BeNllj . submitted). Iterated integrals of 
finite order are Q-functions. 

As in all preceding cases, the proof is constructive and gives an explicit bound for the 
parameters of the Q-function in terms of all relevant integer data, this time including 
the order of the iterated integrals. The bound is given by a tower function (iterated 
exponential) of height 4 and is obviously very excessive. 

4.4. Some open questions. In a somewhat surprising way, all nontrivial examples of 
Q-functions originate from periods of rational forms and their generalizations (iterated 
integrals). It would be very interesting and instructive to have other types of examples, 
not directly related to periods or derived constructions. 

Another question is motivated by Theorem 14. II A quasiunipotent representation p of 
the fundamental group 7Ti(CP m \ S, •) for an algebraic divisor S can be realized by a 
regular flat connection SI (integrable Pfaffian system) with the preassigned monodromy p. 
The system Q is unique modulo a rational gauge transformation (cf. with p.H]) and under 
some conditions (to be determined, depending on the representation and the divisor) is 
defined over Q, producing thus a family of Q-functions. Theorem 14.11 asserts that two 
globally gauge equivalent systems admit rather close (in the relative scale, of course) 
counting functions. The question is whether the counting function can be assigned to the 
pair (ft, S) of the monodromy group and the polar divisor, rather than to a Q-system 
realizing this group. 

Next, we note that the quasiunipotence assumption of Theorem 12.131 can be relaxed: 
in fact, it is sufficient to require that all eigenvalues of the small monodromy operators 
have unit modulus. However, it may well be that such "weakly quasiunipotent" systems 
are necessarily quasiunipotent in the standard (strong) sense (e.g., by some application 
of the Gelfond-Schneider theorem) . 

Finally, the counting problem can be reformulated in geometric terms which admit 
a number of natural generalizations. A regular flat system on CP™ with n-dimensional 
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fibers defines a (singular) foliation & with leaves of dimension m and codimension n 
which, because of the linearity, induces another foliation of dimension m and codimension 
n — 1 on the product <CP m x CP™ -1 . The counting function measures the number of 
isolated intersections between the leaves of this foliation and the (linear or projective, 
respectively) subspaces of the product, which have a complementary (co)dimension n and 
very special position: projection of those subspaces on the base CP m is one-dimensional, 
hence their intersection with the fibers is abnormally large. 

One can try to drop this restriction and ask about the maximal number of iso- 
lated intersections with arbitrary (linear or projective) subspace of the complementary 
(co) dimension. Analytically this means counting isolated solutions of systems of equations 
of Q-functions. This problem seems to be quite challenging and may find applications in 
various counting problems. Besides, one can generalize the settings and consider arbi- 
trary (not necessarily linear) foliations & defined by polynomial data and count isolated 
intersections of their leaves with affine subspaces of complimentary dimension. So far the 
problem is solved only for foliations of dimension 1 [NY99J and of codimension 1 }K91j 
or reducible to the latter case ("Pfaffian chains"). 
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